Activity statistics, avalanche kinetics, and velocity correlations in surface growth 
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We investigate the complex spatio-temporal dynamics in avalanche driven surface growth by means of scaling 
theory. We study local activity statistics, avalanche kinetics, and temporal correlations in the global interface 
velocity, obtaining different scaling relationships among the involved critical exponents depending on how far 
from or close to a critical point the system is. Our scaling arguments are very general and connect local and 
global magnitudes through several scaling relationships. We expect our results to be applicable in a wide range 
of systems exhibiting interface kinetic roughening driven by avalanches of local activity, either critical or not. 
As an example we apply the scaling theory to analyze avalanches and roughening of forced-flow imbibition 
fronts in excellent agreement with phase-field numerical simulations. 
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I. INTRODUCTION 

The dynamics of rough surfaces in systems exhibiting 
avalanches or bursts of activity in response to a slow external 
driving has attracted a great deal of attention in recent years. 
Physical examples include fracture cracks [1], fluid imbibition 
fronts in porous media 

lift 

or the motion of domain-walls 
in magnets [5|]- By means of high resolution fast imaging 
techniques it has recently been possible to characterize exper- 
imentally avalanche dynamics of activity from either the local 
velocity map 01 or the average velocity time series 0, |(| . A 
fundamental problem is to understand the role played by lo- 
cal fluctuations (avalanches) in the building up of large-scale 
space-time interface correlations and scale-invariant roughen- 
ing. 

In this paper we show that local activity statistics, avalanche 
scaling properties, and time correlations (in both, the sur- 
face and the global velocity-fluctuations) are generically con- 
nected through scaling relations that we derive by means of 
simple scaling arguments. Depending on the form of proba- 
bility distribution of the first-return time of activity to a given 
site, we divide avalanche driven surface growth systems into 
two separate classes as follows. We consider that driven sys- 
tems can exhibit either off-critical or critical (fractal) spatio- 
temporal activity, and the scaling relations among the corre- 
sponding exponents are different for each class. On the one 
hand, there are systems where scale-invariant interface growth 
emerges from off-critical spatio-temporal patterns of activity, 
where avalanches have a finite extent and the probability for 
the activity to return to any given place is short tailed. On the 
other hand, there are growing interfaces where local activity 
is critical, showing fractal patterns in space and/or time. Ex- 
amples include systems at the depinning transition as well as 
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interface systems that self-organize into a critical state. This 
leads to critical avalanches of activity with fat tails that lead 
to strong fluctuations of the avalanche size, avalanche extent, 
and first return time of activity. In both cases, however, we 
find that the interface space-time correlations turn out to be 
related to the spatio-temporal activity, leading then to a con- 
nection between local and global quantities via scaling rela- 
tionships that we obtain. Our arguments are very general and 
we expect our results to be widely applicable in systems where 
surface roughness fluctuations are driven by local avalanches. 

As an example of much current interest, we apply the scal- 
ing theory to the problem of forced-flow imbibition fronts, 
which occurs when a viscous fluid advances through a disor- 
dered media displacing a less viscous fluid (typically air) 0], 
In forced-flow imbibition experiments a constant fluid injec- 
tion rate is applied while the spatially averaged velocity of the 
liquid-air interface v is kept constant [4]. This system exhibits 
a natural characteristic length scale £ x ~ (1/v) 1 / 2 0-tOll 
corresponding to the typical extent of interface correlations, 
which can be controlled by the injection rate. Therefore, im- 
bibition offers an excellent example where the theory we shall 
present here can be tested, since the system may be either in 
an off-critical regime (v 3> 0) with avalanches extending over 
short length scales of order £ x , or fronts can be driven into a 
critical regime (V — > 0) by simply tuning the fluid injection 
rate. Interestingly, these regimes can be accesible to experi- 
mentation 14J. 

The paper is organized as follows. In Sec. II we study the 
local aspects of the interface dynamics in terms of activity 
statistics, where we give the exact definitions of critical and 
off-critical activity. We also study there the avalanche kinet- 
ics, where we present the different scaling relations among 
the local quantities describing avalanches. In Sec. Ill we an- 
alyze the global dynamics of the interface by means of the 
multiscaling properties of the surface and the mean interface 
velocity time-correlations. All the different quantities defined 
thoughout the paper are then studied numerically in forced- 
flow imbibition fronts in Sec. IV. The final conclusions are 
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given in Sec. V. 

H. LOCAL DYNAMICS 
A. Activity statistics 

We consider an advancing interface /i(x, t) in d + 1 dimen- 
sions where x is the substrate position and t is time. A site x 
is said to be active at a given time if it is moving. For discrete 
models this corresponds to sites to be updated, h — > h + 1, 
in the next time step. For continuous systems defining single 
site activity is more cumbersome but it can be effectively done 
by considering a site is active whenever the local velocity is 
above some fixed threshold. In order to characterize the lo- 
cal spatio-temporal activity we calculate the first-return time 
probability density, V{(T), for a site x to become active again 
after a period of inactivity T. This probability describes the 
time duration of intervals separating subsequent returns of ac- 
tivity at any given site so that the average number of returns 
M{T) in a time interval of duration T satisfies lfl2ll 

N{T) = T-N{T)( uVi{u)du, (1) 

J ST 

where the temporal resolution can be set to unity, ST = 1, 
without loss of generality. 

In systems far from a critical point, which we denote as 
off-critical systems, activity is expected to be described by a 
general scaling function: 

V { (T)=a(T )f(T/T ), (2) 

such that f(y) —> as y — > oo and a (To) is a normalization 
constant so that f°° duVf(u) = 1 is fulfilled. The average 
return time in this case is simply 

/oo 
uVi(u)du oc T , (3) 

and thus a finite average surface velocity v oc 1/ (T) exists. A 
very common example are systems where the return of activity 
is a Poison process, so that the time interval distribution is 
simply exponential /(T/Tq) ~ exp (—T/Tq). In the case of 
systems exhibiting some degree of correlation between two 
consecutive returns of activity to a given site the distribution 
is more likely to be described by 

P f (T) ^T- T *'exp(-r/r ), (4) 

where t[ represents an effective exponent characterizing the 
power-law regime. This case will also be classified as off- 
critical here since there still exists a finite average surface ve- 
locity with the average return time also given by Eq. ||3}. 

In contrast, under certain circumstances, some systems may 
show patterns of activity recurrence that become critical or 
fractal (i.e. scale-invariant) in time in such a way that the first- 
return time probability density exhibits a power-law asymp- 
totic tail: 

Vf{T) ~ T~ ri for T > 1, (5) 



with exponent Tf > 1 for the distribution to be properly nor- 
malizable. In this case, in the long time limit we obtain that 
the average number of return points is Af(T) ~ T Tt ~ x with 
some critical exponent within the range Tf £ (1,2]. Note that, 
for Tf > 2, although the activity is still described by a power- 
law, the system is actually off-critical due to the existence of 
a finite return time (T), i.e., returns become dense in time, 
and the interface advances at finite velocity (drM ~ const.). 
In contrast, for 1 < Tf < 2 returns rarely occur and activity 
is said to be fractal in time. In this case, the infinite average 
return time 

/oo 
u l ~ Ti du -> oo, (6) 

typical of a critical dynamics. Fractal activity occurs, for in- 
stance, near a pinning/depinning critical transition (v —> 0) 
or in surface growth models that self-organize into the critical 
state. 



B. Avalanche kinetics 

An avalanche is defined as a spatially connected cluster 
of active sites. Avalanches are characterized by the typical 
size s(£) of an event of lateral spatial extent I. The duration 
and lateral extent of an event are related, I ~ i 1 / 2 '™, via the 
avalanche dynamic exponent z av . One expects the average 
avalanche size to scale with the lateral extent up to the cutoff 
length scale (if any): s(£) ~ £ d °- v for £ <c £ x > where d av is the 
avalanche dimension exponent. Conversely, one also expects 
to observe a scaling relation s(t) ~ V between size s and du- 
ration t of avalanches. Since an avalanche of lateral length £ 
leads to a fluctuation of the surface height of the order of the 
local width w(l), d av can be easily related to the local rough- 
ness exponent a\ oc - Indeed, one has s(£) ~ £ d w(£) ~ £ d £ aioc , 
and d av = d+ai oc for a surface in d+1 dimensions. Likewise, 
the avalanche size-duration exponent 7 can then be obtained 
from the identity s ~ t 1 ~ £ d + a ^ and we find 



which connects avalanche dynamics with surface rough- 
ness [2]. 

In the case of off-critical activity, non-critical avalanches 
with a typical spatial extent cutoff £ x are expected so that 
durations are also bounded and scaling only holds for t <C 
i x ~ £x av - The avalanche size can be written in terms of the 
average velocity associated with the event as s ~ £ d Af(t) ~ 
£ d v(£)t, for an avalanche of lateral extent £ and duration t. 
We assume that the typical velocity of the avalanche scales 
with the avalanche size as v(£) ~ £~ b . The identity s ~ 
£ d v{£)t - £ d + a ^ leads to £ z '™~ s - £ ai °° and we obtain the 
scaling relation 

^av = "loc + S, (8) 

relating surface roughness and avalanche exponents in the 
case of off-critical activity (non-critical avalanches). The ex- 
ponent 5 in Eq. (0 gives account of possible inter-event cor- 
relations. For instance, in the case of uncorrelated avalanches, 
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typically, one should expect Gaussian fluctuations of the 
avalanche velocity and (v(£)v(£')) ~ 6 d (£ — £'), giving 
v(£) ~ £~ d / 2 . Avalanches are initiated in places favored given 
the disordered background. If this disorder is spatially uncor- 
rected it follows that avalanche velocity should also be un- 
correlated and 5 = d/2. Correlations in the random medium 
can produce non trivial correlations in the avalanche velocity 
and that possibility is contemplated here. Also, in the case of 
non disordered systems, where movement is driven by some 
external noise, the same argument applies since the noise may 
be either white or correlated. We have assumed a general situ- 
ation by introducing the exponent S that accounts for possible 
velocity correlations in case they exist. 

A different result is obtained for driven surfaces exhibiting 
a fractal profile of activity which leads to critical avalanches 
with no characteristic length or time scales (< x — > oo). 
This is an interesting case because it includes surfaces near 
the pinning/depinning transition, surfaces driven into self- 
organized critical states, and surfaces obeying extremal dy- 
namics, among other surface models in the critical state. As 
discussed above, avalanches in the presence of fractal activity 
do not have a finite velocity since the average number of re- 
turns of activity to a given site is M(t) ~ t Tt ~ x with T{ < 2. 
This means that the scaling relation © does not apply. How- 
ever, the scaling relation between avalanche size and local 
width, s ~ £ d N(t) ~ £ d w(£), is indeed valid also in this case, 
but with Af(t) ~ t 7 * -1 . Thus we find the scaling relation 01311 



Tf 



«loc 



(9) 



which connects the local roughness with the statistics of the 
first-return time of surfaces exhibiting a fractal activity. 



III. GLOBAL DYNAMICS 
A. Multiscaling of the height-height correlations 

Having activity in the form of localized avalanches im- 
plies that growth is highly inhomogeneous in space with many 
time scales involved, which, in turn, may produce multiscal- 
ing of the height-height correlations lfl4ll . As it is shown be- 
low, this can give information about local properties through 
global observables. To this end we define Ah(x,t;to) = 
h(x, t+to) — h(x, to) and investigate the generalized g-height- 
height correlation function in the saturated regime 



C q (t) = (\Ah(x,t;t ) 



-i/g 

Ah(x,t;t )\ q ) 



(10) 



where the over bar is a spatial average and brackets denote 
average over independent realizations. When the exponents 
j3 q depend on q the surface is said to exhibit multiscaling and 
this indicates a highly non trivial probability distribution of 
the height fluctuations. We are interested in obtaining f3 q as a 
function of the roughness and activity exponents. 

In the limit g^oo only the site with the maximum growth, 
A/i max (t) = (ma,x x {Ah(x,t;to)}), contributes and we have 



Coo(t) ~ Ah max (t). Note that Ah max (t) becomes an impor- 
tant quantity since it corresponds to the typical height fluctua- 
tion produced on the surface by avalanche events in a time in- 
terval t. For time intervals t shorter than the average avalanche 
duration cutoff (t <C t x ) we have Aft, max (i) ~ t^°° , where 
the exponent is (3^ = lim^oo (3 q . In contrast, in time 
scales longer than the typical event duration (t ^> t x ) sur- 
face fluctuations are no longer spatially localized — vary many 
avalanches occurring at different parts of the surface overlap 
and contribute to the fluctuation. Therefore surface height 
fluctuations must become spatially homogeneous at long time 
scales {t ^ t x ) and the typical growth A/i max (i) should scale 
with the global surface width W(t) ~ t a ' Za *, where a is the 
global roughness exponent. These two limiting behaviors can 
be put as 



Ah max (t) 



for t < ty 
for t > ty 



(11) 



where the crossover time t x is the average avalanche duration. 
In the case of off-critical activity t x is finite, while it diverges 
for fractal activity. 

For a system of lateral extent L, surface fluctuations over 
time scales t t x involve localized avalanche events of ac- 
tivity and so only a fraction of the surface n mov ~ (£/L) d ~ 
(t 1 / 2av /L) d is moving. By using Eq. ( TTTb . we can roughly 
estimate that a term of the form (Ah) m for m > 1 scales as 



(Ah)™ ~ n mov (t)[Ah maK (t)} m - t m ^ +d '^L' 



(12) 



for fluctuations at time scales f < t x , while it becomes 

(Ah) m - t m(2a+rf)/2z '™, (13) 

for t 3> t x . Now, note that in the temporal regime t <C L z "- V 
of interest here, the main contribution in Eq. ( fTOb comes from 

-1/9 



the term (Ah)i . Therefore we can write C q (t) ~ t^ q , with 



J Ax, + for t « t 

I 2z 



(14) 



that gives the multiscaling exponents /3 q induced in the sur- 
face fluctuations by the localized avalanche dynamics, valid 
for both off-critical (finite t x ) and fractal (unbounded t x ) dis- 
tributions of activity. The exponent for the infinite moment 
can then be obtained by considering that a height fluctuation 
at time scales t <C t x corresponds to an avalanche of lateral 
extent £ and is given by Aft, max ~ s(t)/£ d ~ t~i- d l z ^ . Using 
Eq. O we have (3^ = a loc /z av . 

It is interesting to particularize these results for the case 
of models exhibiting extremal dynamics where the 

avalanche-duration exponent is 7 = 1 (only one site is grow- 
ing at any given time) and we have A/i max (t) ~ s(t)/£ d ~ 
t/£ d and Ax, = 1 — d/z av with t x — > 00 since the model is 
critical, in agreement with existing results for critical depin- 
ning models with extremal dynamics lfl4ll . 

The presence of multiscaling is due to the infinitely many 
time scales involved in the problem that are induced by a 
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highly inhomogeneous surface growth. An illustrative way to 
show this is to note that the time scale defined by C\ (t) ~ t^ 1 
is associated with the avalanche size-duration exponent 7 
since, from Eq. O, Pi = Poo + d/ z av = (a loc + d)/z av = 7 
[cf. Eq. ©]. In the other limit, when q — > 00, we have /Jqo = 
ckioc/^av that corresponds to the time scale of single-site ac- 
tivity. This observation is transparent for a fractal activity dis- 
tribution where A/i max ~ t^°° can be related to the average 
number of returns to a given site, Ah max ~ M(t) ~ f T,_1 . 
This gives Tf = (3^ + 1, which corresponds to the previously 
derived scaling relation This allows us to immediately as- 
sociate the 00- correlation moment (a global observable) with 
the distribution of the return times of activity at any given site 
(a local observable). This is useful for experiments, where 
single site motion, and therefore local observables, may be 
difficult to measure. 



B. Correlations of the global velocity fluctuations 

Recently, there has been much interest (0, 0, S, 03, EH in 
obtaining information about the avalanche kinetics from the 
scaling behavior of the global activity/velocity fluctuations. 
This is particularly important in experiments where the global 
velocity time series is readily available, while it is difficult 
to monitor microscopic avalanches. The scaling relations ob- 
tained in the preceding sections can be used to give a full so- 
lution to this problem in the context of surface growth. 

In the stationary state, to 3> L Zav , we consider the velocity- 
velocity correlation Av(t) — ([v(t + to) — v(to)] 2 ), where 
v(t) is the average instantaneous velocity time series. Equiv- 
alently, we can compute the velocity fluctuation spectrum 
{\v((jj)\ 2 ), where v(w) is the Fourier component at frequency 
uj ~ t . Generally, one expects Av(t) ~ t 2H where H is 
Hurst exponent of the velocity time series, while in Fourier 
space we have {\v(lu)\ 2 ) ~ with v = 2H + 1. Since 

v(t) = dth, the velocity fluctuation spectrum can be written 

as (\v(uj)\ 2 ) = lj 2 (\h(oj)\ 2 ) . We are interested here in the 
case of surfaces exhibiting kinetic surface roughening in d + 1 
dimensions where the structure factor describing the height- 
height correlations scales as 



(\h(k,t)\ 2 ) = fc-( 2Q+d V(^ 1/Zav ), 



where k = |k|, and the scaling function f(u) 
u — ► 0. We have 



leading to the scaling exponent 



Um<|/i(k,u;)| 2 ) 

k— >0 



{l/u 



(2a+d)/z av + l 



2a + d 



1, 



(15) 
vd for 

(16) 
(17) 



which relates the low-frequency behavior of the velocity- 
fluctuation spectrum to the roughness properties of the sur- 
face. This formula generalizes an earlier analytical result lfl9tl . 
v = (d + 4) Jz — 3, valid for the Kardar-Parisi-Zhang equa- 
tion Il20ll . This scaling behavior holds for frequencies above 
~ (l/^) 2av , below which the spectrum becomes flat. 
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FIG. 1. (Color online) (a) First-return time probability density Vi(T) 
calculated at different velocities. It reflects an off-critical activity 
(moderate velocities: v = vo and v — «o/5), characterized by : 

Vt (T) ~ T~ T 'f exp (-T/To) with effective exponents of r} ~ 1.0 
and 1.37, respectively. Fractal activity is found at low velocities: 
v = uo/40 where V s iT) - T~ T f with r f = 1.48 ± 0.02. All 
results shown correspond to a threshold value of cth = 6. (b) Ef- 
fect of the election of a velocity threshold Cth- First-return time 
probability density for the case of v = Vo/40, when the single- 
site activity is critical, calculated by choosing different values of 
the threshold Cth- For cth < 6 the distribution can be fitted to 
V t (T) ~ T' T f exp(-T/r t h) with T th depending on c th . Note 
that Tf — 1.48 is robust for different velocity thresholds cth indicat- 
ing it is a critical exponent. 



IV. CASE STUDY: FLUID IMBIBITION IN RANDOM 
MEDIA 

We now apply the above scaling theory to the problem of 
forced-flow imbibition in disordered media tM- As it has been 
noted in the introduction, the main point in forced-flow imbi- 
bition lays in the existence of a natural characteristic length 
£ x ~ (1/w) 1 / 2 for the typical avalanche extent |0-Qj]], which 
can be controlled by the liquid flow rate. In this way, we have 
a system whose activity can be made progressively fractal as 
the velocity is tuned from moderate to very low values. 

We simulate forced-imbibition in 1 + 1 dimensions in the 
capillary dominated regime by integrating a standard phase- 
field model J3-E1 



-e 2 V 2 0-rKr)], 



(18) 



in a two-dimensional system of lateral size L = 512 with 
e = 1. The quenched random field r](r) > models capillary 
disorder and favors the liquid (wet) phase, forcing the inter- 
face to advance at the expense of the air (dry) phase. The pa- 
rameter M in the above equation is the mobility which we take 
constant at the liquid phase (</> > 0) and zero at the air phase 
(4> < 0). In our numerical model we have used a spatially dis- 
tributed dichotomic quenched noise and the results have been 
averaged over 25 disorder realizations. The interface position 
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FIG. 2. (Color online) Height-height time correlation function 
C q (t) ~ t /3<! , calculated at two typical velocities, shows the presence 
of multiscaling with fa « 0.93, #2 ~ 0.81, /3 3 « 0.73, 04 « 0.70, 
and /Soo w 0.64 at moderate velocities, w = «o (a), and /?i w 0.94, 
/3 2 ~ 0.75, /3 3 = 0.67, /3 4 ~ 0.63, and « 0.53 at low velocities, 
v = v /40 (b). 



FIG. 3. (Color online) Power spectrum of the mean interface veloc- 
ity (|€(w)| 2 ) ~ uj~ v calculated at two typical velocities, obtaining 
v = 1.26(8) at high velocities and v = 1.05(7) at low velocities in 
agreement with the theoretical prediction 1.33 and 1.0 respectively. 



h(x, t) separating the wet and dry phases is computed. A site 
on the interface is active at time t if its local velocity is above 
some fixed threshold, v(x, t) > Ct^v, where cth is some arbi- 
trary constant. 

In Fig [TJa) we plot the distribution of first-return times of 
activity to any given site. We see that for moderate velocities, 
v = vq = 2 x 10~ 3 and v — Urj/5, the activity is an ex- 
ponentially cut-off power-law with an exponent that depends 
on the mean velocity, while it tends to be pure power-law dis- 
tributed for very low velocities, Ti = vq / 40, where we obtain 
Tf = 1.48 ± 0.02 over a range of more than three decades. 
These distributions of activity correspond to the fact that the 
characteristic avalanche size (and duration) diverges as one 
gets closer to the pinned state, 7J — ifTTll . It is worth to men- 
tion here that the observed results do not depend on the choice 
of the arbitrary threshold cth- Figure |TJb) shows that chang- 
ing the value of c t h only introduces an artificial cut-off but 
does not change the exponent, something already observed in 
Ref. JH. 

Let us now study these two different dynamical regimes. In 
the case of moderate velocities (off-critical activity) the sta- 
tionary surface is super-rough [21] and recent numerical work 
shows IH [Till that the global and local roughness exponents 
are, a m 5/4 and a\ oc = 1, respectively. Avalanches in im- 
bibition are driven by the capillary disorder [2], which is a 
random quenched and uncorrelated field, so we should expect 
the Gaussian value 6=1/2. Inserting these values in Eqs. © 
and (8j, we find the dynamic exponent z av = 3/2 and the 
avalanche duration exponent 7 = 4/3 in ci = 1, in excel- 
lent agreement with the existing numerical estimates J2, [HI] 
for imbibition fronts in this range of velocities. Figure [2] 
shows that the surface fluctuations display multiscaling up 
to a given characteristic time scale. In particular, we predict 
Poo = ct\ oc / z av = 2/3, which is in excellent agreement with 
the numerical estimate in Fig. |2ja). However, Eq. ( TBI with 
ot\ oc — 1 and z av = 3/2 only reproduces the trend of the 
numerical estimated exponents but not the exact values, as a 
clear consequence of the strong crossover effects due to the 
finite cut-off time scale t x . We expect our prediction for (3 q 
to become better as t x becomes larger (see below). 

Regarding the velocity fluctuations spectrum we predict a 
decay as ~ (l/w)" with v given by Eq. ( fTTI i up to a cut- 



off frequency uq ~ (l/£x) 2av associated with the satura- 
tion time in a finite-length correlated system. Replacing the 
known numerical estimates J^, [Till for the critical exponents 
a = 5/4 and z av = 3/2, we predict a velocity correlation 
exponent v — 4/3 in good agreement with our simulations 
v = 1.26 ± 0.08 in Fig. [3] as well as the numerical estimate 
in Ref. 0]. 

The scaling behavior changes dramatically when we drive 
the front at very low velocities, where the activity statistics be- 
comes approximatively fractal (see Fig. [T). Although, for an 
infinite size system, a truly fractal activity distribution would 
only appear in the limit v — > 0, we observe that for veloc- 
ities as low as v = Vo/A0 our finite system has a robust 
scale-invariant activity distribution V{ (T) ~ T~ Ti over four 
decades in T. We recall that, according to the arguments dis- 
cussed in the previous sections, the scaling relation ([8]l is no 
longer valid for fractal activity and has to be replaced by Q. 
The scaling exponents a = 3/2, ai oc = 1, z av = 2, and 
7 = 1 in this close-to-pinning regime have been recently ob- 
tained [11] both analytically and numerically by using the 
phase-field model described above. We are therefore able 
to check our scaling relations (O, ( [T4l . and (1171 1. Replacing 
ai oc = 1 and z av = 2 in (O we find the first-return time dis- 
tribution exponent Tf = 3/2, which is in excellent agreement 
with our numerical estimate in Fig.Q] 

Figure |2|b) shows that the system also exhibits multiscal- 
ing in the regime of low velocities. In this regime the system 
is very close to the pinning critical point (lacking character- 
istic lengths or time scales) and our prediction for the mul- 
tiscaling exponents j3 q is expected to be more accurate. In 
fact, inserting a\ oc — 1 and z av = 2 in Eq. ( fl~4b we obtain 
(3 q = (q + l)/2<7, which leads to the set of theoretical ex- 
ponents fa =l,/3 2 = 0.75, p 3 = 0.667, and /3 4 = 0.625 in 
excellent agreement with our simulations [cf. Fig.[2b)]. Also, 
in the limit q — > 00 we expect to have fir*, = a\ oc /z av = 1/2 
to be compared with the numerical result P^ k 0.53. 

Finally, regarding the temporal correlations of the global 
velocity signal in the low-velocity regime we can replace a = 
3/2 and z av = 2 in Eq. ( fl"7l > to obtain v = 1 for the full range 
of frequencies, which is to be compared with our numerical 
estimate v = 1.05 ± 0.07 shown in Fig. [3] 
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V. CONCLUSION 



To summarize, we have presented a scaling theory of sur- 
face roughening in systems driven by avalanches exhibiting 
either off-critical activity distribution, described by a power- 
law with a cut-off function, or critical activity distribution, 
characterized by a pure power-law. Our results connect, on the 
one hand, distributions of local activity with the scaling prop- 
erties of avalanches and surface roughness. On the other hand, 
temporal correlations in both the surface and the global veloc- 
ity signal with local activity and avalanche kinetics turn out 
to be also connected by means of scaling relations. For illus- 
tration, we have applied the theory to the problem of forced- 
fluid imbibition, which can show both off-critical or critical 
distributions of local activity as the the front velocity is tuned 



toward zero, obtaining in both cases excellent agreement be- 
tween theory and numerical results. Our predictions for the 
sigle-site activity statistics, avalanche distribution, roughness, 
and velocity-fluctuation correlations are novel and of interest 
for the experiments in this little explored region in the con- 
text of forced-imbibition near the depinning transition. Our 
arguments are very general and should be of general valid- 
ity in very different contexts where surface roughness fluctu- 
ations are driven by local avalanches of motion, characterized 
by having either a parallel or extremal dynamics. 
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